Sub-wavelength imaging with a left-handed material flat lens 
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We study numerically, by means of the pseudospectral time-domain method, the unique features 
of imaging by a flat lens made of a left-handed metamaterial that possesses the property of negative 
refraction. We confirm the earlier finding that a left-handed flat lens can provide near-perfect 
imaging of a point source and a pair of point sources with clear evidence of the sub-wavelength 
resolution. We illustrate the limitation of the resolution in the time-integrated image due to the 
presence of surface waves. 
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The unique properties of left-handed (LH) materi- 
als 0, i.e. materials with simultaneously negative real 
parts of the dielectric permittivity and the magnetic 
permeability fir (which can also be described by a neg- 
ative index of refraction), allow focusing of electromag- 
netic waves by a flat slab of the material; this effect is in 
sharp contrast to conventional optical lenses with a pos- 
itive refractive index that need to have curved surfaces 
to form an image. Recently, Pendry argued that a 
slab of a lossless LH material with = /^r = ^ 1 should 
behave as a perfect lens enabling one to obtain an ideal 
image of a point source through the reconstitution of the 
evanescent wave components. While recent experiments 
confirmed the main features of negative refraction 0, 0| , 
the question of near-perfect imaging of a flat lens and 
near-field focusing has remained highly controversial 0, 
and is severely constrained because of large dissipation 
and anisotropy in the metamaterial. Nevertheless, sev- 
eral studies showed that nearly-perfect imaging should be 
expected even under realistic conditions when both dis- 
persion and losses of the left-handed material are taken 
into account 0, 0, H H US . 



In this Letter, we re- visit the problem of nearly perfect 
imaging by a flat lens made of a left-handed metamate- 
rial and study numerically, by use of the pseudospec- 
tral time-domain method, imaging by a flat lens made of 
this material. In order to study the amplification of the 
evanescent waves, we compare the wave- vector spectra of 
the field at the image plane in the case with and with- 
out the LH slab. We confirm the earlier finding that a 
left-handed flat lens can provide near-perfect imaging of 
a point source and a pair of point sources with clear 
evidence of sub-wavelength resolution. In addition we 
consider the time-resolved and time-integrated Poynting 
vector at the source frequency which behaves equivalently 
to the time-integrated intensity that may be relevant in 
lithography applications. 

We model an electrical line current source in front of 
a slab of LH material embedded in free space, as shown 
in Fig. ^ The system is translationally invariant in the 
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FIG. 1: Schematic view of the model system. A current 
source is placed in front of a LH slab. In y direction the sys- 
tem models an open domain through the use of reflectionless 
absorbing boundaries (PML), while it is periodic in the x di- 
rection. The solid dot indicates the source location while the 
dashed lines denote the observation planes, i.e. the source and 
image planes. 



z direction and is treated as two-dimensional. The simu- 
lations are performed in TE polarization (H^, Hy, Ez ^ 
0). In the y direction the system is considered open, 
as achieved through reflectionless absorbing boundaries, 
while the x direction is taken as periodic. Because of 
this the simulation essentially uses an infinite slab and 
an array of sources. 

It has been shown that a material with left-handed 
character can only be achieved through dispersion 
Therefore, the LH material is assumed to have lossy 
Drude characteristics in both electric permittivity and 
magnetic permeability 0], given by 
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are the plasma frequencies and 7e, 



2 



are the collision frequencies of the electric and magnetic 
properties, respectively. To simplify and to impedance 
match the slab to free space, we take ujpe — i^pm — and 
7e = 7m = 7- The material parameters are chosen to give 
a refractive index with real part Re(n) = — 1 at frequency 
/o = 15 GHz (cjQ = 27r/o). For this we use ujp = 27r\/2/o. 
The collision frequency is 7 = 27rx4.5 MHz, which results 
in er{i-oo) = /^^(wo) = — 1 + 0.0006*. 

We directly simulate the field propagation based on 
Maxwell's equations, using the pseudospectral time- 
domain method pTl | . In this method all field components 
are sampled at the same spatial location which avoids 
ambiguities of the material interfaces |l2j | inherent in the 
usual finite-difference time-domain method [TsL IT^ , and 
the problems caused by transition layers at the interfaces 
[m . The domain walls in the y direction are covered with 
a uniaxial perfectly-matched layer (PML) 0, 0, 
boundary to simulate an open domain. In the x direc- 
tion the system is periodic. The material constitutive 
relations are implemented using the auxiliary differential 
equation method fl^. 

The current source is turned on slowly using a temporal 
dependence (1— e"*/"^) sin(wot), with a turn-on parameter 
T « 22//0, to reduce its band width and allow the system 
to reach steady-state faster. 

The simulation uses a time step of At = 29.2873 ns and 
a spatial step size of Aa; — Aj^ — Ao/lOO = 0.199862 m; 
Aq is the free-space propagating wavelength at frequency 
/o. The simulation size is 1024 x 128 and is iterated for 
600000 time steps, i.e. 2635 periods, to ensure steady- 
state. 

It has been shown that a slab of LH material with 
parallel sides should focus propagating (]| and evanescent 
[3 waves. Thus, if the slab thickness d is greater than the 
distance a between the object (source) and the front face 
of the slab, then one expects an image to form inside the 
slab, a distance a from the front face, as well as behind 
the slab, a distance d — a from the back face, as indicated 
in Fig. CI 

First, we place a single source a distance a = G.2Ao 
in front of a slab with thickness Ao/3. Snapshots of the 
electric field in the source plane and the two image planes 
are shown in Fig. |2 The snapshots were taken when the 
field was at its maximum in the source plane, which was 
also the time of maximum field in the image planes. The 
field of the source plane is well reproduced in the im- 
age planes. The feature size of the central peak in the 
image planes is well below the free-space limit of Ao/2. 
To illustrate that point we performed a spatial Fourier 
transform of this data and show the transverse spatial 
frequency dependence of the image-plane fields relative 
to the source-plane field in Fig. O For comparison the 
figure also shows the spectrum in the same planes when 
the LH slab is replaced by air. In the case without the 
LH slab, waves with spatial frequency greater than I/Aq, 
corresponding to evanescent waves, are almost entirely 
absent in the image planes. In contrast, when the LH 
slab is present, waves with spatial frequency up to 3.5/Ao 
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FIG. 2: Snapshots of the electric field in the source plane 
(solid), the image plane inside the slab (dotted) and beyond 
the slab (dashed) with a single source. 
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FIG. 3: Transverse wave vector dependence of the fields in 
the image planes relative to the field in the source plane in 
the case of a single source; with LH slab: image location 
inside slab (dotted), image location beyond slab (dashed); 
free-space: image location inside slab (dash-dotted), image 
location beyond slab (solid). 



are transmitted to the image locations. This agrees with 
the smaller than Aq /2 features present in Fig. |21 The 
wave vector cutoff value (fee = 27rAc) seen in Fig. 01 is 
comparable to the value expected according to the pa- 
rameters used in the system At 2.8/Ao a peak occurs 
in the spectrum, similar to the results reported in 20] . 

We also study the imaging of a pair of sources in order 
to characterize the possibility of sub-wavelength resolu- 
tion. The single source is replaced by a pair, separated 
by 0.34Ao. The two sources have the same time depen- 
dence, i.e., are in-phase. In Fig. 0] snapshots of the y 
component of the Poynting vector {Sy) in the image plane 
beyond the slab are shown for two different slab thick- 
nesses while the distance between the sources and the 
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FIG. 4; Snapshot of the y component of the Poynting vector 
in the image plane beyond the slab for a slab thickness of Ao/3 
(solid) and Xq/2 (dashed). 
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FIG. 5: Absolute magnitude of the electric field on a loga- 
rithmic scale at time t = 2635//o with two sources in front of 
the LH slab of thickness Ao/3. The LHM slab is denoted by 
solid lines, the image planes by dashed lines, and the source 
plane by a dash-dotted line. 

slab is keept constant. The snapshots are again taken at 
the same time and show the instance with highest im- 
age peaks. One finds that features with sizes less than 
Ao/2 are reproduced in the images but the quahty of the 
sub- wavelength resolution diminishes with increased slab 
thickness P|. 

With slab thickness Ao/3, Sy in the region between the 
image peaks has opposite sign compared to that of the 
image peaks. This is a signature of the influence of the 
surface waves on the image. 

Figure[51shows a two-dimensional false-color plot of the 
absolute magnitude of the electric field with two sources 
on a logarithmic scale. The geometry and orientation is 
as indicated in Fig. ^ One clearly sees the two sources 
and the strong fields they excite at the surfaces. The 
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FIG. 6: Time dependence of the amplitude of the y compo- 
nent of the Poyting vector (Sy) in the image plane beyond a 
slab of thickness Ao/3 on a linear scale. 



surface waves rapidly diminish within Ao/2 on cither side 
of the source-image axis. One notices that in the image 
planes the source images are resolved in the transverse, 
i.e. X, direction but not in the y direction. Again, we 
see the additional peak between the image peaks and 
extrema of the surface waves to both sides. 

To understand the behavior of the surfave waves we 
studied the transfer function of the system, which is eas- 
ily found by matching the boundary conditions of the 
fields at the slab surfaces. With the parameters of the 
system given above, we found that the appearance of the 
surface waves is dominated by the high wave vector com- 
ponents contained in the source distribution. These wave 
vector components lie beyond the maximum wave vector 
value kc that the slab can image correctly. A source dis- 
tribution with less contribution beyond fee suffers less dis- 
tortion in the image from surface waves. The slab thick- 
ness and the material parameters determine kc which in 
turn influences the occurrance of the surface waves. 

Unlike the cavity effect observed for a slab of finite 
width 21], we found that the distance to the neighboring 
sources (present due to the periodic boundary conditions 
in x) has no noticable influence on the image, as long as 
this distance is at least a few Ac. 

To illustrate the temporal behavior of the image, in 
Fig. we plot Sy in the image plane beyond the slab 
as a function of time for a slab of thickness Ao/3. The 
image peaks are visible as areas of large positive value. 
Between the image peaks, with a time delay of about T/S, 
a negative peak appears. The curve in Fig. ^corresponds 
to a slice at time 0.22r in Fig. El Upon careful inspection 
one can see the peaks of the surface waves moving away 
from the image peaks and then disappear. When the slab 
thickness is increased to Ao/2 (not shown), the region of 
negative Sy has essentially disappeared and become a 
region of relatively large positive Sy. Rather than being 
clearly distinct, the two image peaks begin to merge, as 
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can be seen in Fig. 0] The surface waves present have 
lower kc, i.e., larger Ac, and greater amplitude. They 
also move faster away from the image peaks. 

From Fig. it is evident that the time-integrated 
Poynting vector has an additional significant extremum 
between the two image peaks, which will lead to a dis- 
tortion of the image. In general, this peak due to sur- 
face waves is not distinguishable from a genuinely imaged 
peak in the source plane. This poses an important limi- 
tation on the usable parameter range of the left-handed 
material flat lens. 

In conclusion, we have studied the unique features of 
imaging by a flat lens made of a left-handed metamate- 
rial that possesses the property of negative refraction. By 
employing the pseudospectral time-domain method and 
comparing with the wave propagation in air, we have 
confirmed previous results of near-perfect imaging of a 



point source and a pair of point sources with clear evi- 
dence of sub-wavelength resolution. We found that the 
time-integrated image has an additional limiting factor 
on the resolution of the images of multiple close objects 
due to substantial constructive interference in the image 
plane. We believe that a potential advantage of this kind 
of imaging is its scalability to sub-micrometer scales to 
make possible imaging and nano-photolithography with 
spatial resolution in the tens of nanometers. 
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